Percolation, renormalization, and quantum computing with non-deterministic gates 
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We apply a notion of static renormalization to the preparation of entangled states for quantum computing, 
exploiting ideas from percolation theory. Such a strategy yields a novel way to cope with the randomness of 
non-deterministic quantum gates. This is most relevant in the context of optical architectures, where proba- 
bilistic gates are common, and cold atoms in optical lattices, where hole defects occur. We demonstrate how 
to efficiently construct cluster states without the need for rerouting, thereby avoiding a massive amount of 
conditional dynamics; we furthermore show that except for a single layer of gates during the preparation, all 
subsequent operations can be shifted to the final adapted single qubit measurements. Remarkably, cluster state 
preparation is achieved using essentially the same scaling in resources as if deterministic gates were available. 



PACS numbers: 

In addition to its conceptual interest, the cluster state or 
one-way model of quantum computation UJ appears to yield 
a highly desirable route to quantum computing for a variety 
of technologies [1-7], not least due to the clear cut distinc- 
tion between the creation and consumption of entanglement. 
There are generally two approaches to preparing cluster states: 
In a static approach, one can make use of a physical setting 
where an underlying lattice structure is naturally given, as in 
case of cold atoms in optical lattices. Here, imperfections such 
as hole defects are a challenge for the preparation of a perfect 
cluster. On the other hand, there is a dynamic approach of 
building up large-scale cluster states using probabilistic quan- 
tum gates; this is most promising for architectures based on 
linear optical systems, optical cavities, or optical small non- 
linearities [5-8]. Any such scheme requires dynamics that 
depend on success or failure of the entangling gates. While 
cluster state computation always requires a level of "classi- 
cal" feed-forward - wherein settings of single-qubit measure- 
ment devices need to be switched according to previously ob- 
tained outcomes - all current proposals for building cluster 
states in the dynamic approach rely on larger (by several or- 
ders of magnitude) amounts of the much more daunting "ac- 
tive switching" type of feed-forward. This involves the quan- 
tum systems being routed on demand into different possible 
coherent interactions with other quantum systems, based on 
success or failure of previous gates. This is particularly true 
for the linear optical paradigm [6-10], which is the motivation 
for some of our results, although we stress they are applicable 
to any such probabilistic setting. 

In this work we demonstrate that the two obstacles men- 
tioned above - local defects and active feedforward - can be 
overcome in principle by a single strategy: The appropriate 
use of a static structure together with classical percolation 
ideas. In the dynamic approach, it is possible to dispense with 
all of the active switching, once small initial pieces of clus- 
ter state have been obtained. Given such small clusters, every 
qubit is only involved in one probabilistic two-qubit gate, fol- 
lowed by one single qubit measurement. The principal idea is 
to use the probabilistic gates to combine small pieces of clus- 
ter according to a specially chosen lattice geometry. On the 
percolated lattice lUlll a pattern of single qubit measurements 




FIG. 1: Renormalization procedure: Blocks Ax{k) of the lattice U 
(here shown with overlapping blocks using dashed lines) with cross- 
ing clusters give rise to renormalized sites x £ M. 



can be efficiently determined by an offline classical compu- 
tation, and universal quantum computation is attainable. Re- 
markably, it is possible to achieve this complete removal of 
active feed-forward at essentially no cost: The resources re- 
quired induce at most a sublinear overhead per qubit com- 
pared to the situation of having perfect deterministic gates at 
hand. We will also present strong numerical evidence of a re- 
duction of the overhead to being sub-logarithmic. For linear 
optical settings, we will show how the initial entangled states 
can be as small as 4-qubit cluster states, which have already 
been prepared in down conversion experiments lfl2|] . 

The technique we use to deal with the randomness of the 
cluster states is that of coarse graining of an underlying lat- 
tice U into blocks corresponding to logical qubits, forming a 
renormalized lattice M. Here, vertices comprise the blocks, 
and edges reflect connections between crossing clusters in 
neighboring blocks, see Fig.[T] We want M to be a fully occu- 
pied lattice with asymptotic certainty, and we seek to identify 
the scaling of the resources required to achieve this. 

For concreteness we focus on M = [1,L]^^ for some 
length L, that is, the renormalized lattice is a 2d cubic (square) 
lattice. Of this, a hexagonal sublattice will be used for quan- 
tum computing, the graph states of which constitute universal 
resources. We focus on bond percolation, so a bond is present 
("open") with probability p. The underlying lattice has to be 
chosen such that in some way it is possible to exceed the crit- 
ical probabilities, marking the aiTival of infinite connections 
throughout the lattice 1 1 1], with the initial resources and gates 
at hand. When p already exceeds the critical bond-percolation 
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FIG. 2: (a) Placing 7-qubit clusters at the vertices of a cubic lattice 
and implementing a probabilistic parity check gate (such as a linear 
optical Type-I fusion gate 0]) results in a percolated cluster, (b) For 
quantum computation it suffices to use the 6-qubit graph states 01,01 
depicted (i.e., the complete graph Kq) forming the covering lattice. 

probability pc of a two-dimensional lattice (e.g., 1/2 for the 
square lattice), this natural geometry can be introduced. Then, 
a possible renormalization amounts to simply exploiting verti- 
cal and horizontal paths which by standard results necessarily 
have to cross sufficiently many times. Our goal is to go sig- 
nificantly further and deal also with the situation of small p: 
Techniques to increase p (e.g., replace each bond by multi- 
ple ones in parallel to increase the probability that at least one 
of them exisits) or to decrease p^ (using another lattice with 
higher threshold, notably ones in three or higher dimensions 
to generate a 2d lattice) will be used. This will in general in- 
crease the vertex degree, causing a trade-off between the num- 
ber and the size of the initial cluster state resources. Within 
all solutions which provide efficient scaling of the resource 
number, the main interest will be to minimize their sizes. 

In the most prominent context at hand, namely fusion gates 
^ operating with a success probability of at most 1/2, we 
will overcome this problem by taking U C Z"^, so starting 
from a 3d cubic lattice, for which pc = 0.249. We will iden- 
tify each vertex x G M with a block of size (2fc)'^'^. We can 
meaningfully define an event 'Qlx{k) of x G M being "oc- 
cupied". With this we mean that there exists a crossing open 
cluster within the block, so a connected path on the graph con- 
necting each pair of faces on opposite sides, at least in the first 
and second dimension ifTHl . Moreover, this crossing cluster is 
connected to each of the ones of the blocks associated with 
sites y adjacent to x, which does not arise as naturally as in 
the 2d case. We show the following: 

Renormalized cubic lattices: Let p > p^. Then for any 
/i > 0, the probability Pp{L) of having '^x{k) satisfied for all 
x e M with k = Lf^ fulfills 

lim Pp(L) = 1. (1) 

L — 'oo 

In other words, with a sublinear overhead k = 0{L^^), one 
can create a cubic lattice AI = [1, L] out of U using bond 
percolation. Moreover, this preparation is asymptotically cer- 
tain (in the same sense as in Refs. |§1), despite the underlying 
elements being probabilistic. The value of k specifies to what 
extent we "dilute" the superlattice M compared to U . 

To show the validity of ([T]), we introduce a series of blocks 



of the underlying lattice J7, which, in addition to the blocks of 
M include blocks overlapping with those (see dashed lines in 
Fig.QJ. For any y e [2,2LY'^Mi AJk) = [yik,yik + 2k- 
1] X [ysfc, y2k + 2fc - 1] X [1, 2k] Each vertex x & M 
is identified with y ~ 2x. To show that Pp(2lx(A:)) = 1 (al- 
most certainly) for all x £ M for large L, we make use of 
statements on crossing clusters in cubic lattices, as well as of 
a convenient tool in percolation theory, the FKG inequality: 
Let 5B and £ be two increasing events, i.e., events that "be- 
come more likely" for increasing p. Then the FKG inequality 
states that Pp(QS (1 €) > Pp{T')Pp{€) d. In other words, 
increasing events are positively correlated. 

Let us denote with ^y{k) the event that Ay{k) has a left- 
to-right crossing cluster in the first dimension, i.e., an open 
path having vertices a and b satisfying ai = yik and bi ~ 
yik + 2k ~ 1. Now there exists a constant d > 0, only de- 
pendent on p, such that Pp(5Bj,(fc)) > 1 — cxp(— dfc^) for 
k > 3 ifTlll . We only need to "connect these vertices". The 
blocks Ay{k) and Az{k) are overlapping for dist(?/, z) = 1. 
Now take two sites y e [2, 2iy — 1] x [2, 2L], with y2 even, 
and z with zi = yi + 1, and Z2 = 2/2- The events ^y{k) 
and ^z{k) are increasing events, and therefore, we can use 
the FGK inequality: intuitively, if in Ay{k) there is already 
a crossing cluster, then this crossing cluster is already half 
way through Az{k), and hence renders a crossing cluster 
there more likely. Consider the overlap between two adjacent 
blocks, By{k) = Ay{k){^Az (fc). We can define the following 
event: For p S [0, Pc], we define Dy [k) as the event that never 
occurs, forp G {p^, 1] it is the event of having at most a single 
left-to-right crossing cluster in By{k). This is an increasing 
event lUSll . Hence, the probability of having simultaneously a 
left-to-right crossing cluster in Ay{k), one in Az{k), and ex- 
actly one in By (fc) can be estimated using the FKG inequality. 
There exist constants c, a > 0, only dependent on p such that 
the probability of having the event 'Z)y{k) satisfies lITsIl 

Ppi^yik)) > 1 - {2kfaexp{-ck). (2) 

So, using again the FKG inequality, one finds that the proba- 
bility of 2;j,(fc) of having two crossing clusters in Ay{k) and 
^(i;i+i,y2)(^) which are actually connected as Pp{€y{k)) > 
(1 — exp(— o?fc^))^ (1 — {2k)^a exp(— cfc)). This procedure can 
be iterated, using FKG in each step. To find connections in 
the other direction, we re-use the argument on having at most 
a single crossing cluster, but now using [1, 4A:|jn the third di- 
rection to be able to apply the results of Ref. lUSll . This gives 
an overall probability of having ^x{k) for each x € M of 

PpiL,k)^Pp{Ux^xm > (l-exp(-dfc2))2L^-Lx 

((1 - (2fc)6acxp(-c/c))2(l - (4/c)6acxp(-c2fc)))^^-^"^' . 

We bound this expression from below with the slowest in- 
creasing term, i.e., there exists an integer ko such that 
Pp{L,k) > (1 - {2k)^aexp{-ck))^^'' for all k > 
ko. Let us set k = for e > {). Then, 

ImiL^oo Pp{L,k{L)) = 1 using that for any e, / > 0, 
we have that lim„^oo (l — en'^^ cxp(— /n'^/^))" = 1. This 
means that by using a sublinear overhead, we arrive at an 
asymptotically certain preparation of the renormalized lattice. 
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This gives rise to an overall resource requirement of 
0{Uf X = 0(L2+3e-) j.qubit states to build a fully con- 
nected cluster state that (almost certainly) consists of L x L 
blocks, and requires no rerouting. As long as p > p^, this scal- 
ing will hold. This should be compared to the 0{L^) qubits 
we would require if we had perfect deterministic gates with 
which to build the cluster 

To utilize the renormalized blocks some classical computa- 
tion is needed, and we need to ascertain that it is efficient in the 
system size. One first has to find the crossing clusters in each 
block, e.g., by the Hoshen-Kopelman-algorithm ||16i l {0{k^) 
steps, 0{k^) memory). Then, using a series (0(L^)) of 
breadth-first-searches {0{k^)) equips us with loop-free con- 
nections between chosen "mid-points" of all sets of neighbor- 
ing blocks. For conceptual simplicity it suffices to only iden- 
tify 3-way "T-junctions" (instead of crosses), and use these to 
build a hexagonal lattice, which can be converted to a square 
lattice easily ifisll . The unused qubits in the blocks can be 
eliminated using az measurements, all unwanted qubits on the 
connections between the T-junctions can be cut out with 
and ay measurements. Thus, the amount of classical compu- 
tation needed shows the same scaling behavior as the number 
of required resources lfl4ll . 

At this stage we have used 7-qubit clusters on a cubic lat- 
tice, see Fig.|2] We now turn to various methods for reducing 
the size of this initial resource. The first one is quite gen- 
eral, and will apply to any lattice. We see from Fig. |2]that 
a qubit is left on each successfully formed edge. One inter- 
esting observation is that this qubit may be measured out, re- 
laxing the requirement of photon number resolving detectors 
to dichotomic detectors. However, one might also use this to 
construct the covering lattice [ITT] of the original lattice, by 
connecting these sites with all perimeter sites from the neigh- 
boring stars, and removing the stars' central qubits (Fig.|2b)). 
From percolation theory it is known ifTHl that the critical bond 
percolation probability of a lattice equals the critical site per- 
colation probability of the covering lattice (for which a site is 
"open" with a certain probability p). Thus by using 6-qubit 
clusters (with the connectivity of the complete graph Kq as 
shown) the covering lattice can be built by fusion of neighbor- 
ing corner qubits. These percolation processes are equivalent 
for our purposes, because a path between two arms of one star 
in the original lattice exists iff the fusion processes involving 
these two arms were successful, and a path between two cor- 
ner qubits in the covering lattice exists iff the fusion attempts 
on the equivalent two qubits were successful. 

A quite different method (somewhat more specific to linear 
optics) can further reduce the size of the initial states on the 
cubic lattice to stars with 3 arms where the central qubit is 
redundantly encoded by judiciously fusing the two "central" 
qubits of each of the 5-stars, while simultaneously applying 
the Type-I fusion operations on the bonding qubits (see Fig. 
|3]l. A Type-I and a Type-II fusion are applied to the centers. 
On success of one fusion gate, the centers merged into a sin- 
gle redundantly encoded qubit and subsequent application of 
another fusion gate will succeed and only reduce the level of 
encoding. If the first fusion fails the second one may still suc- 
ceed with p ~ 1/2. As the order does not matter, both gates 




FIG. 3: A pair of 5-qubit states (star with 3 arms, central qubit re- 
dundantly encoded) can be used to create a single 7-qubit GHZ state 
with p = 3/4. A Type-I and a Type-II fusion are applied to the cen- 
ters. On success of one fusion gate, the centers merged into a single 
redundantly encoded qubit and subsequent application of another fu- 
sion gate will succeed and only reduce the level of encoding. If the 
first fusion fails the second one may still succeed with p = 1 /2. As 
the order does not matter, both gates may be applied simultaneously, 
without any need for rerouting. The bond fusion processes are not 
influenced, only the site probability is determined by this process. 



may be applied simultaneously, without any need for rerout- 
ing, so site-preparation succeeds with probability 3/4. If the 
central ("site") fusion fails, the bond fusions can still be at- 
tempted as usual. The single qubits resulting from the failure 
are in the state |+)®^, and fusion gates involving them will 
succeed or fail with probability 1/2. Hence, the site and edge 
generation processes are independent and do not require active 
switching, allowing for application of the mixed percolation 
model fllll. 

A more general approach to decrease the size of the ini- 
tial resources is the following 11711 : Instead of using the cu- 
bic lattice, we switch to a different 3d lattice, e.g., the one 
with the lowest vertex degree, namely the diamond lattice 
which has vertex degree 4, and a bond percolation threshold 
of pc = 0.389. While percolating on the diamond lattice di- 
rectly would require 5-qubit star clusters, by percolating on 
the covering lattice (as explained above), the pyrochlore lat- 
tice, we even further reduce the initial resources required to 
4-qubit GHZ states, which lies below the resource size dic- 
tated by a naive ansatz with a square lattice! These tetrahedra 
consist of triangles and are thus not two-colorable. However, 
it can readily be shown that the resulting graph states can still 
be reduced to universal cluster states. 

As less is known analytically about percolation for the di- 
amond lattice, we have turned to a numerical verification that 
this lattice suffices for our purposes. In fact, we find that the 
resource scaling appears slightly more favorable than the up- 
per bound proven above for the cubic lattice. Cubic blocks 
of the diamond lattice of size k^^ have been simulated and 
arranged in two dimensions as described above, then used as 
renormahzed lattice. These sites are occupied iff there exist 
crossing clusters connecting the four faces. Bonds between 
neighboring sites exist iff the crossing clusters of the corre- 
sponding blocks are connected through the common face. De- 
pending on k and the probabilities of a site and an edge being 
open, the probability P{L) of building up the whole renormal- 
ized lattice of size L x L without any missing sites or bonds 
is obtained. By requiring a fixed threshold P, the scaling of 
the block size k{L) that is needed to lie above this threshold is 
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FIG. 4: Dependence of the diamond lattice block size k^'' on the size 
L of the renormalized square lattice for for three different sets of site- 
and bond probabilities (psitcPbond). The overall success probability 
threshold P{L) was chosen to be 1/2. 10''' blocks of each size were 
created; each lattice site was randomly populated 10'^ times. 

found. The results are summarized in Fig.|4] which suggests a 
scaling of k{L) ~ o(log(i)) for each set of parameters, thus 
a scaling of L^o{log (L)) of 4-qubit cluster states to build a 
lattice of size L x L with a success probability of at least P. 

Although the primary concern of this work is to show how 
to effectively remove active feed-forward and to deal with 
probabilistic gates or configurations using percolation tools, 
we briefly address the role of further imperfections in such 
a setting. To start with, losses, the occurrence of which is 
known without the need for destructively measuring them 
(e.g., atomic systems with probabilistic entangling gates fl), 
can be accounted for by measuring out the qubits around fail- 
ure sites in the cr^-basis. This is effectively a new percolating 
model, albeit with correlated probabilities. Numerical inves- 
tigations show that a similar reasoning as above is expected 
to hold, demonstrating that losses of 10% can easily be ac- 
counted for. More importantly, concerning losses like pho- 
ton loss, which are only detected by destructive qubit mea- 
surements, this scheme is no different than others, i.e., stan- 
dard techniques can make it loss tolerant, although, needless 
to say, with significantly more effort. Our approach read- 
ily suggests two strategies to cope with such errors: On the 
one hand, by fixing the block size, one would fix the effective 
(non-heralded) site occupation probability. Then, schemes for 
fully-fletched fault-tolerant one-way computation can be used, 
once above the respective fault-tolerant threshold liTsll . On the 
other hand, to suppress loss rates, specifically photon loss, it 
is legitimate to consider initial encodings like tree-structures 
ifisll . They in turn can be grown probabilistically as well, still 
without any need for active switching. For example, a crude 



estimate is the following: To correct for these losses occuring 
with a rate of 10%, trees with branching parameters (6, 7, 7, 1) 
can be used to suppress the loss to an effective rate of 10"^. 
Together with blocks of size 6^'^ elementary diamond cells 
(which fail to be crossed with a probability of 4.76 • 10^^), 
the overall site loss rate on the renormalized lattice lies below 
the 3 -10-30] limit 10]. 

In this work we have introduced a method based on perco- 
lation phenomena of building cluster states with probabilis- 
tic entangling gates. The scheme dramatically reduces the 
amount of coherent feed-forward required; specifically there 
are no rerouting steps needed, once one starts from appropri- 
ate building blocks which can be as small as 4-qubit states. 
We provided a proof that to prepare an L x L cluster state, 
asymptotically with certainty, even with this very restricted 
set of tools, a scaling in the number of resources of 0{L^^^) 
for any e > can be achieved. Numerical simulations have 
been carried out, suggesting an even better resource consump- 
tion of _L^o(log'3(i)), which should be compared to in the 
case of perfect deterministic gates. 

We emphasize that while we have focussed on a rigorous 
proof of the utility of one type of percolation model, there 
are many obvious ways in which our work can be extended: 
different lattices, site percolation, dynamically growing per- 
colated clusters and finding optimal, loss corrected paths "on 
the fly", and much more. We expect therefore our ideas to 
be applicable to a wide range of architectures where proba- 
bilistic quantum gates originate, e.g., from exploiting small 
non-linearities as in Ref. IjStl, or linear optics jLTj]. In the con- 
text of ultracold atoms in optical lattices, cluster states may 
be prepared by exploiting cold collisions One could then 
think of universal computational resources when starting with 
a Mott state exhibiting hole defects, such that the filling is not 
exactly a single atom per site, giving rise to a site percola- 
tion variant of our argument. It would also be interesting to 
see whether the new freedom of measurement-based schemes 
for quantum computing beyond the one-way computer ll20ll 
gives rise to further improvements concerning resource re- 
quirements. The presented ideas open up a new way to deal 
with randomness of probabilistic gates in quantum computing. 
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